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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  NO.  1678 

THEORETICAL  WAVE  DRAG  AND  LIFT  OF  THIN 
SUPERSONIC  RING  AIRFOILS 
By  Harold  Mirels 


SUMMARY 

An  approximate  linearized  solution  is  presented  for  the  wave 
drag  and  lift  of  an  airfoil  generated  by  rotating  a  thin  supersonic 
profile  about  an  axis  parallel,  or  nearly  parallel,  to  its  chord. 

The  aerodynamic  coefficients  are  obtained  from  a  surface  distribu¬ 
tion  of  sources,  of  strengths  proportional  to  the  local  airfoil 
slopes,  about  a  cylinder  whose  radii’s  and  chord  equal  those  of  the 
original  ring  airfoil.  This  source  distribution  satisfies  the 
boundary  conditions  when  the  part  of  the  wing  within  the  forward 
Mach  cone  from  a  point  on  the  wing  surface  departs  only  slightly 
from  a  plane.  The  solution  is  therefore  accurate  for  ring  airfoils 
having  chords  that  are  small  in  comparison  to  the  radius  of  rotation. 

The  lift  coefficient  of  thin  supersonic  ring  airfoils,  based 
on  the  airfoil-surface  area,  is  one-half  the  Ackeret  value  for  a 
two-dimensional  wing  of  infinite  span.  The  drag  coefficient  is 
equal  to  the  sum  of  the  Ackeret  value  for  the  given  profile  (with 
the  ring  airfoil  at  zero  angle  of  attack)  and  the  induced  drag 
coefficient.  These  coefficients  are  probably  within  5  percent  of 
the  correct  linearized  values  for  ring  airfoils  whose  chord-radius 
parameter  (chord  divided  by  the  product  of  the  radius  nnfl  the 
cotangent  of  the  Mach  angle)  is  within  the  range  from  0  to  0.20. 


INTRODUCTION 

The  linearized  solution  for  the  aerodynamic  coefficients  of  a 
thin  supersonic  airfoil  whose  chord  elements  form  a  flat  surface  is 
currently  in  an  advanced  stage.  Conical  superposition  has  been 
extensively  used  to  determine  performance  when  the  plan  form  is 
bounded  by  straight-line  segments.  Surface  source  distributions 
have  been  used  to  develop  analytical  and  numerical  methods  of  solu¬ 
tion  (references  1  and  2)  for  an  arbitrarily  shaped  plan-form 
boundary.  However,  airfoils  whose  chord  surface  is  not  flat  have 
been  comparatively  neglected. 
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Wave  drag  lift  are  evaluated  herein  for  the  supersonic 
ring  airfoil  generated,  by  rotating  a  thin  sharp-edged  profile 
section  about  an  axis  parallel,  or  nearly  parallel,  to  its  chord. 
The  chord-radius  ratio  of  the  ring  airfoil  is  assumed  to  be  small. 
A  surface  source  distribution,  as  originally  proposed  in  refer¬ 
ence  3  for  essentially  flat  airfoils,  is  used  as  the  basis  for  the 
analysis.  The  limiting  chord-radius  ratios  for  which  this  source 
distribution  gives  a  valid  linearized  solution  are  estimated  by  ^ 
comparing  the  results  for  the  external  surface  of  the  ring  ej1*011 
to  those  obtained  by  the  numerical  method  of  reference  4.  The 
ring  airfoil  in  subsonic  flight  is  considered  in  reference  5. 

The  investigation,  conducted  at  the  NACA  Cleveland  laboratory, 
•was  completed  during  December,  1947 . 


SYMBOLS 


°D 

Cf 


The  following  symbols  are  used  in  this  report: 


drag  coefficient  based  on  wing  surface,  /  -r — - — 

IJpU2  2*rc/ 

skin-friction  drag  coefficient  based  on  wing-surface  area, 
/friction  drag's 
\  |  PE2  2^c  ) 

lift  coefficient  based  on  wing-surface  area,  / - — - ■) 

l^ptr  2*rc/ 


pressure  coefficient, 
wing  chord 


incremental  pressure 

,  ipB2  / 


D  drag 

L  lift 

M  free- stream  Mach  number 

q  local  source  strength  per  unit  area 

r  mean  radius  of  airfoil 
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s  area  of  Integration 

t  maximum  profile  thickness 

t/c  thickness  ratio 

U  free-stream  velocity,  taken  in  positive  i  direction 
v  perturbation  velocity  normal  to  source  surface 


r 


Cartesian  coordinate  system 


i 

r 

6,  co 


>  cylindrical  coordinate  system 


a  angle  of  attack,  radians 

ctQ  angle  between  chord  of  profile  section  and  axis  of  ring 
airfoil,  radians 

0  cotangent  of  Mach  angle,  (Vm2-1) 

A  local  slope  of  wing  surface  with  respect  to  direction  of 
free-stream  velocity 

p  free-stream  density 

a  local  slope  of  wing  surface  with  respect  to  chord 

cp  perturbation-velocity  potential 

Subscripts: 

e  external  surface  of  ring  airfoil 

i  internal  surface  of  ring  airfoil 
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METHOD  OF  ANALYSIS 

The  linearized  equation  for  the  perturbation-velocity  potential 
of  an  irrotational,  compressible  fluid  is  given  by 


(i  - 


aftp  b2* 

Sr2  dy2  dz‘ 


A  general  solution  of  equation  (l)  can  be  obtained  from  a  surface 
distribution  of  sources.  The  potential  at  a  point  in  the  flow  field 
is  the  sun  of  the  contributions  of  the  elemental  sources  in  the 
forward  Ifech  cone  from  that  point.  If  the  source  surface  is  in  the 
x,y  plane,  the  potential  at  a  point  in  the  x,y  plane  is  given  by 


vhere  q  is  the  function  defining  the  source-strength  distribution. 

The  coordinate  system  and  the  limits  of  integration  are  illustrated 
in  figure  l(a) . 

The  induced  normal  perturbation  velocity  at  a  point  on,  or  an 
infinitesimal  distance  from,  a  plane  containing  sources  is  simply 
a  function  of  the  local  source  strength  and  equals  q*.  (See 
reference  3.)  Thus,  if  a  planar  distribution  of  sources  is  to 
represent  the  flow  about  the  top  or  bottom  surface  of  a  thin  super¬ 
sonic  airfoil  whose  chord  surface  is  flat,  the  boundary  conditions 
are  satisfied  by  setting 

q  »  ]£  =  22  (3) 

a  Jt 

where  X  =  v/u  is  the  local  slope  of  the  wing  surface  in  the  free- 
stream  direction.  If  the  sources  are  assumed  to  be  in  the  x,y  plane, 
the  equation  for  the  velocity  potential  becomes 


(4) 
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The  following  sign  convention  for  A  will  he  used:  When  the  surface 
normal  projecting  into  the  stream  has  a  component  in  the  positive 
x  direction,  the  slope  is  considered  negative;  when  the  surface 
normal  has  a  component  in  the  negative  x  direction,  the  slope  is  con¬ 
sidered  positive.  For  a  flat  plate  at  angle  of  attack  a,  the  slopes 
for  the  upper  and  lower  surfaces  are  constant  and  have  the  values 
A=  -a  and  A=  a,  respectively. 

The  potential  flow  about  a  thin-ring  airfoil  can  he  found  from 
a  source  distribution  about  a  circular  cylinder  whose  axis  is  par¬ 
allel  to  the  free-stream  direction  and  whose  radius  and  chord  equal 
those  of  the  original  airfoil;  however,  the  function  defining  the 
source- strength  distribution  remains  to  be  determined.  In  the 
limiting  case  of  a  thin-ring  airfoil  having  a  small  chord-radius 
ratio,  the  wing  section  in  the  forward  Mach  cone  from  a  point  on  the 
wing  departs  only  slightly  from  a  plane.  From  the  previous  dis¬ 
cussion,  a  source  distribution  defined  by  the  local  wing  slope 
q  =  Ml/jt  adequately  satisfies  the  boundary  conditions.  This  rep¬ 
resentation  is  directly  analogous  to  the  solution  of  an  essentially 
flat  wing  in  that  allowance  is  made  for  variations  of  wing  slope, 
but  the  departure  of  the  wing  surface  from  a  plane  is  considered 
negligible.  The  velocity  potential  for  points  on  a  ring  airfoil  of 
small  chord,  expressed  in  cylindrical  coordinates,  is  then 


where  rdw,  cor,  and  0r  have  replaced  dq,  tj,  and  y  of  equa¬ 
tion  (4).  The  coordinate  system  is  illustrated  in  figure  1(b). 

The  normal  perturbation  velocity  at  a  point  on  a  three- 
dimensional  source  surface  is  not  a  simple  function  of  the  local 
source  strength.  Equation  (3),  and  consequently  equation  (5), 
therefore  do  not  satisfy  the  boundary  conditions  when  the  chord- 
radius  ratio  of  a  ring  airfoil  is  sufficiently  large  that  the 
deviation  from  a  plane  of  the  chord  surface  in  the  forward  Mach 
cone  is  no  longer  negligible. 


WAYS  DRAG  AND  LIFT  OF  CYLIimRICAL-RING  AIRFOIL 

A  cylindrical-ring  airfoil  is  a  ring  airfoil  of  zero  profile 
thickness,  camber,  and  flare  (<xq  =  0).  The  velocity  potential 
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of  such  an  airfoil  may  be  found  if  the  function  defining  X  at  all 
points  is  known.  3For  small  angles  of  attack  and  a  small  chord, 
this  function  is  shown  in  appendix  A  to  be 

XQ  =  -  a  sin  a)  (6) 

Ai  =  a  sin  to  (6a) 

The  expression  for  the  velocity  potential  of  the  external  surface 
of  a  cylindrical-ring  airfoil  results  from  the  substitution  of 
equation  (6)  into  equation  (5): 


It  is  desirable  to  solve  equation  (7)  directly  for  dcpe/c*x.  In 
appendix  B,  t'his  solution  shows  that 


The  pressure  coefficient  for  a  point  on  the  external  surface,  based 
on  the  linearized  Bernoulli  equation,  is  then 
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The  pressure  coefficient  for  the  internal  surface  is  the  negative 
of  that  for  the  external  surface. 

The  lift  and  drag  coefficients,  referred  to  the  airfoil- surface 
area,  may  be  obtained  from  the  respective  summations  of  the  effec¬ 
tive  lift  and  drag  forces  acting  on  each  element  of  area: 


The  results  for  the  lift  and  drag  coefficients  show  that 


The  series  in  equations  (9),  (12),  and  (13)  converge  very 
rapidly  for  values  of  x/(pr)  and  c/(p?)  that  are  sufficiently  small 
to  justify  the  use  of  a  source  distribution  defined  by  the  local 
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wing  slopes.  No  appreciable  error  will  be  introduced  by  assuming 
these  series  to  equal  1.  The  pressure  coefficient  then  becomes 
simply  the  Ackeret  value  corresponding  to  the  local  slope: 

C-n  _  =  I  (-a  sin  8)  (14) 

P>e  3 

This  result  indicates  that  the  spanwlse  variation  of  slope  in  the 
forward  Mach  cone,  due  to  chord-surface  curvature,  may  be  neglected 
when  calculating  the  pressure  coefficient  at  points  on  an  airfoil 
by  means  of  a  source  distribution  defined  by  q  =  Mj/jt.  Therefore 
equation  (14)  could  have  been  obtained  directly  from  equation  (7) 
by  using  the  mean  value  sin  0  instead  of  sin  co. 

Similarly,  the  lift  and  drag  coefficients  become 

0L-f  f15) 

and 

°D  -  T  (16) 

Both  coefficients  are  one-half  the  Ackeret  values  for  a  flat-plate 
airfoil  of  infinite  span. 


WAVE  BRAS  AND  LIFT  OF  RING  AIRFOIL  OF  GENERAL  PROFILE 

In  the  general  case,  a  ring-airfoil  profile  section  has  camber 
and  thickness  and  its  chord  makes  a  finite  angle  with  the  axis  of 
the  airfoil.  The  local  airfoil  slopes  are  then  given  by  the 
expressions 


Xe  =  Oe  -  (oq  +  a  sin  8)  (17) 

^i  =  Of  +  (oq  +  a  sin  0)  (17a) 

For  a  small  chord,  the  pressure  coefficient  has  been  shown  to  equal 
approximately  the  Ackeret  value  corresponding  to  the  local  slope. 

cp,e  =  |[°e  -  (otg  +  a,  sin  8)]  (18) 

=  2  [pi  +  (oq  +  a  sin  8)J  (18a) 
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The  lift  coefficient  obtained  by  substituting  these  pressure 
coefficients  in  equation  (10)  equals  that  of  the  cylindrical-ring 
airfoil.  Thus  the  lift  coefficient  of  a  ring  airfoil  is  independent 
of  its  profile  slope. 


The  drag  coefficient  is 


r  c 


(oe2  +  a±2)  dx 


(19) 


The  first  term  of  equation  (19)  defines  the  induced  drag  due  to 
lift  and  is  one-half  the  Ackeret  value  for  induced  drag  coefficient. 
The  last  two  terms  represent  the  drag  of  the  ring  airfoil  at  zero 
angle  of  attack  and  equal  the  drag  coefficient  of  an  Ackeret  airfoil 
(having  the  given  profile  section)  at  angle  of  attack  Oq. 


LIFT-DRAG  RATIOS  OF  RING  AIRFOIL 

A  calculation  can  be  made  to  illustrate  the  order  of  magnitude 
of  lift-drag  ratios  obtainable  from  a  ring-type  airfoil.  For  a 
symmetrical  diamond  profile,  no  flare,  and  a  skin-friction  drag 
coefficient  Cf ,  the  lift-drag  ratio  is  given  by 


(20) 


For  small  values  of  t/c  and  Cf,  the  L/D  performance  approaches 
that  of  an  Ackeret  airfoil.  With  increasing  t/c,  the  L/D 
approaches  one-half  the  value  for  the  corresponding  Ackeret  airfoil. 
Equation  (20)  can  be  maximized  to  give  the  angle  of  attack  at  which 
(L/D)max  occurs  for  a  given  t/c  and  pCf.  The  result  shows 


and  the  corresponding  value  of  (L/D)maI  is 


(21) 
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Equation  (22)  is  plotted  in  figure  2.  The  values  of  airfoil  lift 
coefficient  at  which  occurs  are  cross-plotted  in  this 

figure.  As  shown  by  equations  (15),  (21),  and  (22),  this  lift 
coefficient  is  the  reciprocal  of  (L/D)^ia-r. 

Little  information  is  available  on  the  value  of  friction  drag 
coefficients  at  supersonic  speeds.  The  mean  value  Cp  =  0.006  is 
sometimes  used  for  thin  supersonic  airfoils.  If  C^,  is  assumed 
independent  of  Mach  number,  figure  2  indicates  that  (L/D)^^ 
decreases  with  increased  flight  speed,  but  this  effeot  becomes 
less  pronounced  with  larger  values  of  t/c. 

Because  the  flare  was  considered  zero  and  a  diamond  profile 
has  the  lowest  possible  wave  drag  for  a  given  thickness  ratio,  the 
values  for  (l/EO^x  presented  in  figure  2  are  the  largest  values 
theoretically  obtainable  from  a  ring-type  airfoil. 

VALUES  OF  c/(pr)  FOE  WHICH  AERODYMMIC 
COEFFICIENTS  AEE  ACCURATE 

The  values  of  chord-radius  ratio  beyond  which  equations  (15) 
aryl  (19)  are  no  longer  valid  can  be  estimated  by  comparing  the 
external  lift  and  drag  coefficients  as  indicated  by  these  equations 
to  the  values  found  from  the  numerical  method  of  reference  4,  which 
is  an  extension  of  the  solution  for  slender,  pointed-nose  bodies  of 
revolution  presented  in  references  6  and  7.  A  line  distribution  of 
sources  and  doublets  is  placed  along  the  axis  of  an  open-nosed  body 
of  revolution,  and  the  strength  distribution  is  so  adjusted  that 
the  flow  and  cross-flow  components  of  the  free-stream  velocity 
follow  the  external  contour  of  the  body.  The  resulting  values  for 
external  wave  drag  and  lift  may  be  considered  the  correct  linearized 
solution.  The  internal  flow,  however,  is  not  properly  described. 

By  use  of  the  line-doublet  distribution  of  reference  4,  the 
lift  coefficient  for  the  external  surface  of  a  cylindrical-ring 
airfoil,  expressed  in  the  fora  PC^  e/a,  was  found  to  be  solely  a 
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function  of  the  airfoil  parameter  c/(fir).  The  relation  is  shown 
in  figure  3.  The  decrease  in  lift  with  increase  in  chord  indicates 
that  the  downstream  sections  of  the  cylindrical-ring  airfoil  have 
progressively  lower  external  lift.  Equation  (15),  "based  on  the 
surface  source  distribution,  yields  a  constant  value  for  PCj^e/a, 
but  the  discrepancy  between  the  two  methods  is  less  than  5  percent 
for  values  of  c/(pr)  up  to  0.20. 

Computations  for  the  external  lift  and  drag  of  finite-thickness 
ring  airfoils  indicate  that  the  effect  of  thickness  ratio  on  the 
discrepancy  between  the  two  methods  is  small  for  the  magnitudes  of 
thickness  ratio  permitted  by  the  linearized  theory.  These  discrep¬ 
ancies  depend  primarily  on  the  value  of  c/(|3r) .  This  result  is  to 
be  expected  ,as  chord-plane  curvature,  and  not  thickness  ratio,  is 
the  source  of  the  difference.  If  the  discrepancy  for  the  internal 
surface  is  assumed  to  be  of  the  same  order  as  that  for  the  external 
surface,  the  generalization  may  be  made  that  the  coefficients 
expressed  in  equations  (15)  and  (19)  are  correct  to  within  5  per¬ 
cent  for  values  cf  c/((3r)  up  to  0.20.  The  degree  of  error  may  be 
estimated  by  the  plot  of  figure  3. 


SUMMARY  OF  RESULTS 

An  approximate  linearized  solution  for  the  aerodynamic  coef¬ 
ficients  of  ring  airfoils  having  small  chord-radius  ratios  indicates 
that  the  lift  coefficient,  based  on  the  wing-surface  area,  is  inde¬ 
pendent  of  the  profile  shape  and  equals  one-half  the  Ackeret  value 
for  a  two-dimensional  wing  of  infinite  span.  The  drag  coefficient 
equals  the  Ackeret  value  for  the  given  profile  (with  the  ring  air¬ 
foil  at  zero  angle  of  attack)  plus  the  induced  drag  coefficient. 
These  coefficients  are  probably  within  5  percent  of  the  correct 
linearized  solution  for  values  of  chord-radius  ratio  from  0  to  0.20. 


Flight  Propulsion  Research  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  April  30,  1948. 
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APPENDIX  A 


FUNCTION  DEFINING  LOCAL  SLOPES  ON  CYLINDRICAL  BING  AIRFOIL 

A  cylindrical  ring  airfoil  of  small  chord  and  at  a  small  angle 
of  attach  is  assumed. 


for  which 

Ay  vertical  displacement  between  corresponding  points  on  leading 
and  trailing  edges 

Ar  radial  distance  between  corresponding  points  on  leading  and 
trailing  edges 

Then 

Av  =  c  sin  a 
Ar  =  Av  sin  go 

=  c  sin  a  sin  go 


/ 


But 


c 


=  sin  a  sin  co 


QftQ 
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Inasmuch  as  a  is  small 

A  =  a  sin  co 

and  with  the  proper  sign  convention 

A  e  =  -a  sin  co 

and 


Xi  b  a  sin  co 
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APPENDIX  B 

t 


SOLUTION  POE  FROM  VELOCITY- POTENTIAL  EQUATION 

bx 


From  equation  (7) 


Because  sin  co  is  independent  of  t,  the  order  of  integration 
of  equation  (Bl)  can  he  reversed  and  the  equation  integrated  with 
respect  to  | .  The  resulting  expression  for  <Pe  can  then  be  dif¬ 
ferentiated  with  respect  to  x  to  yield  a  line  integral  equation 
for  b%/bx.  Reference  1  shows  that  the  result  is  equivalent  to  a 
line  integration  along  the  leading  edge  1=0  of  the  airfoil. 


4 


QCQ 
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Let 


e  -  co 

x/(8r) 


Then 

r6+x/(pr) 


6-x/(pr) 


cos  co  sin 


1  fe  -  col 

|_x/(pr)J 


dco  =  -  — 


X 

Pr 


J 


-1 


cos  [9 


(-3) 


ZX^  .  -1  - 

—  ’  sin  x  zdz 


=  -  —  P  ( so s  6  cos  ~  +  sin  9  sin  sin-'*-  zdz 

Pr  Jl  \  3r  Pr/ 


(B4) 


And  by  sine  and  cosine  series  expansions 


=  -  JL  cos  6 

r1 

1  _  _1  ( *x\2  +  JL  (l 

sx\4 

-  ]  “  •  •  • 

pr 

'i 

2i  \  pr/  41  \| 

irj 

sin--*-  zdz 


x  .  .  f"1  /  zx  V  1  /zx\3 

"  p§  8ln  e  Jx  +  *  ‘  ‘J 


sin--*-  zdz  (B5) 


But 


f  zn  sin-1 
J1 


zdz  =  0 


for  n  equal  to  a  positive  even  integer  or  zero,  and 
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Figure  2.  -  Maximum  lift -drag  ratios  obtainable  from  ring 
airfoil  having  symmetrical  diamond  profile  and  no  flare* 
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